Search for: All records

Changes in global mean sea level (GMSL) during the late Cenozoic remain uncertain. We use a reconstruction of changes in δ¹⁸O of seawater to reconstruct GMSL since 4.5 million years ago (Ma) that accounts for temperature-driven changes in the δ¹⁸O of global ice sheets. Between 4.5 and 3 Ma, sea level highstands remained up to 20 m above present whereas the first lowstands below present suggest onset of Northern Hemisphere glaciation at 4 Ma. Intensification of global glaciation occurred from 3 Ma to 2.5 Ma, culminating in lowstands similar to the Last Glacial Maximum lowstand at 21,000 years ago and that reoccurred throughout much of the Pleistocene. We attribute the middle Pleistocene transition in ice sheet variability (1.2 Ma to 0.62 Ma) to modulation of 41-thousand-year (kyr) obliquity forcing by an increase in ~100-kyr CO₂variability. 

Abstract. We use a recent reconstruction of global mean sea surface temperature change relative to preindustrial (ΔGMSST) over the last 4.5 Myr together with independent proxy-based reconstructions of bottom water (ΔBWT) or deep-ocean (ΔDOT) temperatures to infer changes in mean ocean temperature (ΔMOT). Three independent lines of evidence show that the ratio of ΔMOT / ΔGMSST​​​​​​​, which is a measure of ocean heat storage efficiency (HSE), increased from ∼ 0.5 to ∼ 1 during the Middle Pleistocene Transition (MPT, 1.5–0.9 Ma), indicating an increase in ocean heat uptake (OHU) at this time. The first line of evidence comes from global climate models; the second from proxy-based reconstructions of ΔBWT, ΔMOT, and ΔGMSST; and the third from decomposing a global mean benthic δ18O stack (δ18Ob) into its temperature (δ18OT) and seawater (δ18Osw) components. Regarding the latter, we also find that further corrections in benthic δ18O, probably due to some combination of a long-term diagenetic overprint and to the carbonate ion effect, are necessary to explain reconstructed Pliocene sea-level highstands inferred from δ18Osw. We develop a simple conceptual model that invokes an increase in OHU and HSE during the MPT in response to changes in deep-ocean circulation driven largely by surface forcing of the Southern Ocean. Our model accounts for heat uptake and temperature in the non-polar upper ocean (0–2000 m) that is mainly due to wind-driven ventilation, while changes in the deeper ocean (> 2000 m) in both polar and non-polar waters occur due to high-latitude deepwater formation. We propose that deepwater formation was substantially reduced prior to the MPT, effectively decreasing HSE. We attribute these changes in deepwater formation across the MPT to long-term cooling which caused a change starting ∼ 1.5 Ma from a highly stratified Southern Ocean due to warm SSTs and reduced sea-ice extent to a Southern Ocean which, due to colder SSTs and increased sea-ice extent, had a greater vertical exchange of water masses. 

Abstract. To understand the erosivity of the eastern portion of the Laurentide Ice Sheet and the isotopic characteristics of the sediment it transported, we sampled buried sand from deglacial features (eskers and deltas) across eastern Canada (n = 10), a landscape repeatedly covered by the Quebec-Labrador Ice Dome. We measured concentrations of 10Be and 26Al in quartz isolated from the sediment and, after correcting for sub-surface cosmic-ray exposure after Holocene deglaciation, used these results to determine nuclide concentrations at the time the ice sheet deposited the sediment. To determine what percentage of sediment moving through streams and rivers currently draining the field area was derived from incision of thick glacial deposits as opposed to surface erosion, we used 10Be and 26Al as tracers by collecting and analyzing modern river sand sourced from Holocene-exposed landscapes (n = 11). We find that all ten deglacial sediment samples contain measurable concentrations of 10Be and 26Al equivalent on average to several thousand years of surface exposure – after correction, based on sampling depth, for Holocene nuclide production after deposition. Error-weighted averages (1 standard deviation errors) of measured 26Al/10Be ratios for both corrected deglacial (6.1 ± 1.2) and modern sediment samples (6.6 ± 0.5) are slightly lower than the production ratio at high latitudes (7.3 ± 0.3) implying burial and preferential decay of 26Al, the shorter-lived nuclide. However, five deglacial samples collected closer to the center of the former Quebec-Labrador Ice Dome have much lower corrected 26Al/10Be ratios (5.2 ± 0.8) than five samples collected closer to the former ice margins (7.0 ± 0.7). Modern river sand contains on average about 1.75 times the concentration of both nuclides compared to deglacial sediment corrected for Holocene exposure. The ubiquitous presence of 10Be and 26Al in eastern Quebec deglacial sediment is consistent with many older-than-expected exposure ages, reported here and by others, for bedrock outcrops and boulders once covered by the Quebec-Labrador Ice Dome. Together, these data suggest that glacial erosion and sediment transport in eastern Canada were insufficient to remove material containing cosmogenic nuclides produced during prior interglacial periods both from at least some bedrock outcrops and from all glacially transported sediment we sampled. Near the center of the Quebec-Labrador Ice Dome, ratios of 26Al/10Be are below those characteristic of surface production at high latitude. This suggests burial of the glacially transported sediment for at least many hundreds of thousands of years and the possibility that ice at the center of the Quebec-Labrador Ice Dome survived many interglacials when more distal ice melted away. 

Tropical glaciers have retreated over recent decades, but whether the magnitude of this retreat exceeds the bounds of Holocene fluctuations is unclear. We measured cosmogenic beryllium-10 and carbon-14 concentrations in recently exposed bedrock at the margin of four glaciers spanning the tropical Andes to reconstruct their past extents relative to today. Nuclide concentrations are near zero in almost all samples, suggesting that these locations were never exposed during the Holocene. Our data imply that many glaciers in the tropics are probably now smaller than they have been in at least 11,700 years, making the tropics the first large region where this milestone has been documented. 

High-resolution records from past interglacial climates help constrain future responses to global warming, yet are rare. This dataset contains seasonally-resolved climate records from subarctic-Canada using micron-scale measurements of oxygen isotopes (δ18O) in speleothems with apparent annual growth bands from three interglacial periods – Marine Isotope Stages 11 (409-376 ka), 9 (336-305 ka) and MIS 5e (123-118 ka). Our study highlights the potential for high-latitude speleothems to yield detailed isotopic records of Northern Hemisphere interglacial climates beyond the reach of Greenland ice cores and offers a framework for interpreting them. Table S1 contains the Uranium-Thorium dates for six speleothems, or more specifically, flowstones, from a cave in Northwest Territories (NWT), Canada. It also contains constructed age models for each sample. Then, we applied a two-tiered methodological approach to reconstruct past subarctic climate. First, we produce an ultra-high-resolution δ18O record that, although not continuous, spans thousands of years for portions of these interglacials. This record was created using Secondary Ion Mass Spectrometry (SIMS) to measure δ18O approximately every 35-micrometer (µm) down each sample’s growth axis. This data is shown in Table S2. Second, we used Confocal Laser Fluorescence Microscopy (CLFM) to identify several fluorescent annual bands in each speleothem, which we then targeted for additional SIMS measurements. This data is shown in Table S3. Though these subarctic speleothems are small in size (most are less than 10 centimeter (cm) in length), the application of both CLFM and SIMS on these samples demonstrate their potential for providing ultra-high-resolution records of high-latitude Northern Hemisphere terrestrial climate outside of Greenland and provide insights into interpretive frameworks for future cold-region speleothem δ18O records. 

While there are no ice sheets in the Northern Hemisphere outside of Greenland today, it is uncertain whether this was also the case during most other Quaternary interglacials. We show, using in situ cosmogenic nuclides in ice-rafted debris, that the Laurentide Ice Sheet was likely more persistent during Quaternary interglacials than often thought. Low 26Al/10Be ratios (indicative of burial of the source area) in marine core sediment suggest sediment source areas experienced only brief (on the order of thousands of years) and/or infrequent ice-free interglacials over the past million years. These results imply that complete Lauren- tide deglaciation may have only occurred when climate forcings reached levels comparable to those of the early Holocene, making our current interglacial unusual relative to others of the mid-to-late Pleistocene. 

Abstract While there are no ice sheets in the Northern Hemisphere outside of Greenland today, it is uncertain whether this was also the case during most other Quaternary interglacials. We show, using in situ cosmogenic nuclides in ice-rafted debris, that the Laurentide Ice Sheet was likely more persistent during Quaternary interglacials than often thought. Low 26Al/10Be ratios (indicative of burial of the source area) in marine core sediment suggest sediment source areas experienced only brief (on the order of thousands of years) and/or infrequent ice-free interglacials over the past million years. These results imply that complete Laurentide deglaciation may have only occurred when climate forcings reached levels comparable to those of the early Holocene, making our current interglacial unusual relative to others of the mid-to-late Pleistocene. 

Accurate reconstruction of Laurentide Ice Sheet volume changes following the Last Glacial Maximum is critical for understanding ice sheet contribution to sea-level rise, the resulting influence of meltwater on oceanic circulation, and the spatial and temporal patterns of deglaciation. Here, we provide empirical constraints on Laurentide Ice Sheet thinning during the last deglaciation by measuring in situ cosmogenic 10Be in 81 samples collected along vertical transects of nine mountains in the northeastern United States. In conjunction with 107 exposure age samples over five vertical transects from previous studies, we reconstruct ice sheet thinning history. At peripheral sites (within 200 km of the terminal moraine), we find evidence for ∼600 m of thinning between 19.5 ka and 17.5 ka, which is coincident with the slow initial margin retreat indicated by varve records. At locations >400 km north of the terminal moraine, exposure ages above and below 1200 m a.s.l. exhibit different patterns. Ages above this elevation are variable and older, while lower elevation ages are indistinguishable over 800−1000 m elevation ranges, a pattern that suggests a subglacial thermal boundary at ∼1200 m a.s.l. separating erosive, warm-based ice below and polythermal, minimally erosive ice above. Low-elevation ages from up-ice mountains are between 15 ka and 13 ka, which suggests rapid thinning of ∼1000 m coincident with Bølling-Allerød warming. These rates of rapid paleo-ice thinning are comparable to those of other vertical exposure age transects around the world and may have been faster than modern basin-wide thinning rates in Antarctica and Greenland, which suggests that the southeastern Laurentide Ice Sheet was highly sensitive to a warming climate. 

Abstract. There is unambiguous evidence that glaciers have retreated from their 19th century positions, but it is less clear how far glaciers have retreated relative to their long-term Holocene fluctuations. Glaciers in western North America are thought to have advanced from minimum positions in the Early Holocene to maximum positions in the Late Holocene. We assess when four North American glaciers, located between 38–60∘ N, were larger or smaller than their modern (2018–2020 CE) positions during the Holocene. We measured 26 paired cosmogenic in situ 14C and 10Be concentrations in recently exposed proglacial bedrock and applied a Monte Carlo forward model to reconstruct plausible bedrock exposure–burial histories. We find that these glaciers advanced past their modern positions thousands of years apart in the Holocene: a glacier in the Juneau Icefield (BC, Canada) at ∼2 ka, Kokanee Glacier (BC, Canada) at ∼6 ka, and Mammoth Glacier (WY, USA) at ∼1 ka; the fourth glacier, Conness Glacier (CA, USA), was likely larger than its modern position for the duration of the Holocene until present. The disparate Holocene exposure–burial histories are at odds with expectations of similar glacier histories given the presumed shared climate forcings of decreasing Northern Hemisphere summer insolation through the Holocene followed by global greenhouse gas forcing in the industrial era. We hypothesize that the range in histories is the result of unequal amounts of modern retreat relative to each glacier's Holocene maximum position, rather than asynchronous Holocene advance histories. We explore the influence of glacier hypsometry and response time on glacier retreat in the industrial era as a potential cause of the non-uniform burial durations. We also report mean abrasion rates at three of the four glaciers: Juneau Icefield Glacier (0.3±0.3 mm yr−1), Kokanee Glacier (0.04±0.03 mm yr−1), and Mammoth Glacier (0.2±0.2 mm yr−1). 

Much of our understanding of Cenozoic climate is based on the record of δ¹⁸O measured in benthic foraminifera. However, this measurement reflects a combined signal of global temperature and sea level, thus preventing a clear understanding of the interactions and feedbacks of the climate system in causing global temperature change. Our new reconstruction of temperature change over the past 4.5 million years includes two phases of long-term cooling, with the second phase of accelerated cooling during the Middle Pleistocene Transition (1.5 to 0.9 million years ago) being accompanied by a transition from dominant 41,000-year low-amplitude periodicity to dominant 100,000-year high-amplitude periodicity. Changes in the rates of long-term cooling and variability are consistent with changes in the carbon cycle driven initially by geologic processes, followed by additional changes in the Southern Ocean carbon cycle.